Experimental multiparticle entanglement dynamics induced by decoherence 
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Multiparticle entanglement leads to richer correlations than two-particle entanglement and gives 
rise to striking contradictions with local realism [1], inequivalent classes of entanglement [2], and 
applications such as one-way or topological quantum computing [3, 4]. When exposed to decohering 
or dissipative environments, multiparticle entanglement yields subtle dynamical features and access 
to new classes of states and applications. Here, using a string of trapped ions, we experimentally 
characterize the dynamics of entanglement of a multiparticle state under the influence of decoher- 
ence. By embedding an entangled state of four qubits in a decohering environment (via spontaneous 
decay), we observe a rich dynamics crossing distinctive domains: Bell-inequality violation, entangle- 
ment superactivation, bound entanglement, and full separability. We also develop new theoretical 
tools for characterizing entanglement in quantum states. Our techniques to control the environment 
can be used to enable novel quantum-computation, state-engineering, and simulation paradigms 
based on dissipation and decoherence [5-7]. 
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When exposed to an environment, bipartite entan- 
glement already shows subtle dynamical features, e.g., 
finite- time disentanglement [8, 9]. In a multipartite set- 
ting, decoherence and dissipation enable novel quantum 
applications [5-7], induce even more interesting dynam- 
ics due to the different classes of states, and can decrease 
the number of particles genuinely entangled -an effect 
recently observed [10]. However, decoherence can also 
influence many other state properties useful for quantum 
information processing, such as distillability and the en- 
tanglement between subsystems. A state is distillable if, 
using local operations and classical communication, one 
can extract the maximally entangled states required by 
several quantum communication protocols such as dense 
coding [11] and teleportation [12]. 

An environment can also drive a multiparticle en- 
tangled and distillable state into the undistillable but 
still entangled domain [13, 14], called bound entangled 
(BE) [15]. This class of states is expected to appear in 
many body systems [16], and despite of being undistil- 
lable, it is useful for entanglement superactivation [17], 
quantum secret sharing [18], or remote information con- 
centration [19]. BE states have also given insights into 
classical information theory: a classical analog of bound 
entanglement, called bound information, exists [20]. Re- 
cently, a BE state was simulated with photons [21], as 
well as in a variant called pseudo-bound entanglement 
using nuclear magnetic resonance (NMR) [22]. 

Here, we report the experimentally-observed dynamics 
of entanglement and distillability in the neighborhood 
of a BE state under a partially decohering environment. 
Entanglement and distillability of a multiparty system 
are defined with respect to the state bipartitions, or ab- 
stract splits into two subsystems. Our work focuses on 
a four-party system which can be bi-partitioned in two 



ways, either in pairs, or 2:2, and as a single party plus 
the rest, or 1:3. In this case, we call the state of the four 
particles 2:2 {1:3) separable if every 2:2 (1:3) bipartition 
can be written as a mixture of bipartite states IV^/e), 



and a, /3, /i, and v denote the particles. Otherwise, if 
every 2:2 (1:3) bipartition cannot be written as above, 
we call the state 2:2 {1:3) entangled. Regarding distil- 
lability, a state is 2:2 {1:3) distillable if, for every 2:2 
(1:3) bipartition, a Bell pair can be distilled and each 
element of the pair belongs to one subsystem. An even 
stronger distillability property is entanglement superac- 
tivation, which in the case of a four-particle state, it en- 
ables five parties sharing two copies of the state to distill 
entanglement between the two parties holding a single 
particle [17]. 

Our study starts by preparing a 2:2- and l:3-entangled 
state, violating a CHSH-type Bell inequality [18], and 
capable of entanglement superactivation. As we apply 
a tunable decohering environment, the state stops vio- 
lating the Bell inequality. Then, only within a region 
of further decoherence the entanglement superactivation 
protocol is successful, while the 2:2 and 1:3 entangle- 
ment is preserved with even more decoherence. Increas- 
ing the decoherence eventually eliminates the entangle- 
ment in both bipartitions at different but finite times. 
This finite-time disentanglement behaviour is also some- 
times called environment-induced sudden death of en- 
tanglement [8, 9]. Since the 2:2 entanglement disappears 
before the 1:3, we realize a domain, which can be called 
BE [13]. Similarly, a recent theoretical study showed that 
a four-qubit GHZ state (Greenberger-Horne-Zeilinger [1]) 
can decay into a BE state by becoming 2:2 undistillable 



while slightly 1:3 entangled [13]. Further decoherence 
eventually makes the state fully separable, long before 
the single-particle coherence would asymptotically dis- 
appear. 

Our experiment proceeds in three stages: state prepa- 
ration, exposure to tunable decoherence, and state char- 
acterization. The goal of the first stage is to gener- 
ate a state which decays into the domain of BE states 
when exposed to a partially decohering mechanism. We 
chose to prepare an initial state close to the Smolin state 
ps [23], a known four-qubit BE state usually written as 
a mixture of the Beh states ^^ = (|00) ± |ll))/\/2 and 
^± = (|01)±|10))/V2, 

PS = [|^+^+)(-| + |^-^-)(-| + |^+^+)(-| + |^-^-)(-|]/4, 

(1) 
where |0) and |1) are the qubit basis states, and we use 
the notation |x)('| = lx)(xl- The Smolin state can be 
prepared by first noticing that it is a mixture of four 
GHZ-like states 

PS = [(lllll) + |0000))(-| + (11100) + |0011))(-| + 

(11010) + |0101))(-| + (11001) + |0110))(-|]/8 . (2) 

This state can be reached by applying a single-step GHZ- 
entangling operation to the mixture 

p = [|1111){-| + |1100){-| + |1010)(-| + |1001)(-|] /4, (3) 

where the operation takes a state of the form 1x1X2X3X4) 
into (IX1X2X3X4) -|- |xiX2X3X4))/\/2 and Xi denotes the 
complement of state Xi of ion i in the computational ba- 
sis. The mixture described by equation (3), in turn, can 
be generated by completely decohering a state in which 
three out of four particles are entangled, 

1^) = (lllll) + |1100) + |1010) + llOOl)) /2 , (4) 

with the same mechanism as used in the second stage 
of our study. Finally, another GHZ-entangling operation 
and an NMR-like refocussing technique applied to the 
state 1 1111) generates the state in equation (4). 

In the second stage, the intended rich dynamics was 
achieved by increasingly decohering the initial state, as 
described below (see also Supplementary Information). 
We characterized the state's entanglement and distilla- 
bility in the last stage. A single criterion, the Peres- 
Horodecki separability criterion [24] can prove undistill- 
ability while its extension into a measure, known as nega- 
tivity [25], can quantify entanglement. According to this 
criterion, if a state is separable then its partial trans- 
position has no negative eigenvalues (it has a positive 
partial transpose, PPT). On the other hand, it has been 
shown that PPT states are undistillable [15]. Therefore, 
entangled but undistillable states, or BE states, can be 
detected by verifying that every 1:3 bipartition has a neg- 
ative partial transpose (entanglement) [26], while every 



2:2 bipartition has a PPT (undistihability) [15]. To deter- 
mine the state's undistihability properties, we performed 
a complete tomographic reconstruction. Full knowledge 
of the state enabled us to also check further separability 
and distillability properties. Especially, we designed a 
novel algorithm to prove separability of the states, which 
is a stronger statement than undistihability (see Supple- 
mentary Information). 

Our work was performed on a system of four ^^Ca+ 
ions confined to a string by a linear Paul trap with ax- 
ial (radial) vibrational frequencies of approximately 1.2 
MHz (4.4 MHz). Each ion hosts a qubit on the elec- 
tronic Zeeman levels 1)5/2 (^^ = ~l/2), encoding |0), and 
Sii2{Tn = —1/2), encoding |1), determined by a mag- 
netic field of ^ 4 G. The ion string was optically pumped 
to the starting quantum state |1111) after being Doppler 
cooled and sideband cooled to the ground state of the 
axial center-of-mass (COM) mode [27]. The state of the 
qubits can be manipulated via (i) collective unitary op- 
erations U{0^(j)) — exp(— ^l^^), with S<^ — ^^ ^ 



Zl/e=l' 



1^) 



cos(0)cri^ ' + sin(0)cry , and Gn ' a Pauli spin 
operator acting on the /cth ion, (ii) single-qubit light- 
shift operations Z^^\0) = ex.p l—i^az ), and (iii) a 
GHZ-entangling operation known as M0lmer-S0rensen 
gate [28, 29], MS{0,(j)) = exp (-ilSl). We can prepare 

four-qubit GHZ states with a fidelity of 96% and perform 
collective unitaries and light-shift operations at a fidelity 
of 99%. These imperfections determine the proximity of 
our prepared initial state to the Smolin state. The full 
experimental sequence is shown in Fig. 1. 

The partially decohering mechanism indicated in Fig- 
ure 1 was implemented in the four steps shown in Fig- 
ure 2: (i) hiding the population in |0) by a full coher- 
ent transfer into Si/2{m = 1/2); (ii) transfer of the 
population in |1) into the superposition ^/Y^^\l) + 
y^|I)5/2(^ = —5/2)); (iii) quenching of the population 
in Dr^i2{m = —5/2) into P^/2ij^ = ~3/2) by exposure 
to 854-nm radiation, so that it spontaneously decays to 
|1); and finally (iv) restoring the hidden population into 
|0). In this way, a fraction 7 of the population in |1) irre- 
versibly loses phase coherence with |0) by tracing over the 
emitted photon. In this case we call this basis-dependent 
partial loss of coherence decoherence in the |0);|1) basis. 
In our experiment, we decohere the states in the |0) ± |1) 
(|0) ± z|l)) basis by applying the collective unitary ro- 
tations [/(7r/2,7r/2) ([/(7r/2,0)) prior to the above deco- 
herence, as shown in Fig. 1. The complete decohering 
step in the preparation of the intermediate state in equa- 
tion (3) was performed in the computational basis with 
7 = 1. 

The dynamics of entanglement was explored by vary- 
ing the amount of decoherence 7 to which the initial state 
was exposed (see Fig. 1). After being partially decohered, 
the density matrices of the prepared states p(7) were to- 



mographically reconstructed (see Fig. 3). Error analysis 
was performed via Monte Carlo (MC) simulations over 
the raw data outcomes of the state tomography. The 
amount of entanglement and signature of undistillability 
of the measured states as a function of decoherence 7 
are shown in Figure 4; the explicit values for the most 
representative states are quoted in Table 1. Figure 4 
also indicates other properties of the states determined 
independently of the plotted data (full details in Supple- 
mentary Information). 

The measured initial state (7 = 0) is highly entangled 
in the 1:3 bipartitions (7Vi:3 ^ 0) and slightly entangled 
in the 2:2 (A^2:2 > 0); in addition, it violates a CHSH-type 
Bell inequality and is capable of entanglement superac- 
tivation. The properties of the state already change at 
7 = 0.06, when the state no longer violates the tested 
Bell inequality. The entanglement superactivation pro- 
tocol is successful in the domain of states from 7 = 
up to 7 = 0.12. We show strong evidence in the Supple- 
mentary Information that all measured states from 7 = 
to 7 = 0.18 are biseparable. This means that, although 
they are entangled with respect to any fixed 1:3 and 2:2 
bipartition, they can be written as a mixture of separa- 
ble states, which are separable with respect to different 
bipartitions. 

The passage into bound entanglement occurs at 
7 ^ 0.21. While the measured state at 7 = 0.24 is 2:2 
separable and 1:3 entangled, the bound entanglement is 
arguable because a fraction of the MC samples revealed 
2:2 entanglement, thus indicating insufficient statistics. 
By 7 = 0.32, the state is now bona fide BE, when also 
all MC samples are 1:3 entangled, 2:2 undistillable, and 
even 2:2 separable. The change into full undistillability 
is heralded by the state measured at 7 = 0.47 because 
all eigenvalues of the partial transpose were positive for 
every bipartition, as shown in Table I and Figure 4. How- 
ever, while the measured state is fully separable, known 
methods failed to prove the separability of the MC sam- 
ples. By 7 = 0.60, we achieve full separability in the 
measured data and all MC samples. 

We also found a BE state by decohering the initial 
state in only the |0) ± |1) basis, represented by the state 
at y = 0.43 (see Table I). We thus controllably realized a 
passage into bound entanglement in a simple decohering 
environment, with statistically significant entanglement 
in the 2:2 partitions and positivity of the eigenvalues of 
the 1:3 partially-transposed states. 

In conclusion, we experimentally explored the dynam- 
ics of multiparticle entanglement, separability, and dis- 
tillability under a tunable decohering mechanism. The 
influence of the environement naturally created a bound 
entangled state. Our investigation on the dynamics of 
multiparticle entanglement can be extended to observe 
bound entanglement on other states such as decaying 
GHZ states [13] or thermal states of spin models [16]. In 
addition, recent quantum-computing, state-engineering. 



and simulation paradigms driven by dissipative or deco- 
hering environments [5-7] can benefit from the environ- 
ment engineering techniques here demonstrated. 
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TABLE I: Negativity and smallest eigenvalue of the partial transpose of representative states. The parties in the bipartitions 
are labelled as A, B, C, and D. The state with 7^ = 0.43 was prepared by decohering in a single step, see text. Uncertainties 
in parentheses indicate one standard deviation, calculated from propagated statistics in the raw state identification events. 

AB:CD ACiBD ADiBC AiBCD B:ACD CiACD DiABC 



7 = 


A^2:2 = 


0.033(5) 


0.041(5) 


0.044(5) 


A^l:3 


7 = 0.32 


min[eig(p^2:2)] ^ 


0.020(2) 


0.019(2) 


0.022(2) 


A^l:3 


7 = 0.47 


min[eig(p^2:2)] ^ 


0.028(3) 


0.029(3) 


0.031(2) 


min[eig(p^i^3)] 



0.715(8) 0.715(8) 
0.035(7) 0.032(8) 
0.019(3) 0.020(3) 



0.715(8) 0.716(8) 
0.038(8) 0.045(7) 
0.019(3) 0.018(3) 



:0.43 min[eig(p^^-^)] = 0.013(2) 0.011(2) 0.013(2) 



A^i:; 



0.038(7) 0.039(7) 0.042(7) 0.045(7) 



( 


3-out-of-4 
sntangling gate 


full rii_i7 

df£°[- gate 




partial decoherence 

7 7 


P(7) 

t 


ll>- 
ll>- 
ll)- 
ll>- 


CO 


-Z(7r)- 


o" 

^■^ 
^ 

^ 
§ 




0^ 

/ 





0^ " 






f 


0^ 


0" 








Q. 
CD 05 

l| 













































Initial state preparation 



Environment 



FIG. 1: Experimental sequence for studying the dynamics of 
multiparticle entanglement and distillability induced by de- 
coherence. The operations indicated in the sequence are col- 
lective unitary transformations [/(^, 0), light-shift gates ^(^), 
and entangling gates MS{0^ (j)). 
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FIG. 2: Zeeman-split ^°Ca^ levels for the implementation of 
tunable decoherence via entanglement with a spontaneously- 
decaying photon. Partial decoherence is realized by simulta- 
neously performing on all ions the steps: (i) hiding of |0), (ii) 
partial transfer of |1), (iii) quenching and decay of ^/7|1), and 
finally (iv) restoring |0). 



|p(7)l for 7 = 



|p(7)l for 7 = 0.32 



|p(7)l for 7 = 0.60 




FIG. 3: Density matrices (absolute value) of states which are 2:2 and 1:3 entangled (7 = 0), bound entangled (7 = 0.32), and 
fully separable (7 = 0.60). The components of the Smolin state (Eq. 2) are highlighted with distinctive colors. The density 
matrices of all measured states are shown in the Supplementary Information. 
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FIG. 4: Negativity and smallest eigenvalue of the partial transpose for each 2:2 and 1:3 bipartition of the measured states as 
a function of decoherence. A positive smallest eigenvalue of the partial transpose, min[eig(/>^)] > 0, reveals undistillability. 
Bipartitions data are slightly offset horizontally for clarity, but all visible groups correspond to the same amount of decoherence 
indicated by the tick marks. Error bars indicate zb 1 standard deviation, calculated from propagated statistics in the raw state 
identification events. The solid lines were calculated by decohering the prepared initial state with a 0.05 offset in 7 (due to 
imperfections in the decoherence implementation) . The properties shown in bold were determined by tests independent of the 
plotted data (see Supplementary Information). 
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I. PROVING SEPARABILITY OF QUANTUM STATES 

A. A direct algorithm to prove separability for a bipartite 
state 

Before explaining our novel algorithm for proving sep- 
arability of a quantum state, let us first define our no- 
tations and definitions. We consider a bipartite N x M 
system with Hilbert space H = Ha ^ ^b- Any matrix g 
acting on 1-L which is Hermitian [g = ^"f"), has no negative 
eigenvalues {g > 0), and is normalized {Tr{g) = 1) is a 
valid density matrix of some quantum state. 



By definition, a state is separable, if it can be written 
as a convex combination of product states. 



Q = ^Pk\ak){ak\ ^ \hk){^k 



(1) 



where the pk are non-negative {pk ^ 0) and normalized 
C^kPk — 1)5 i^ other words, they form a probability 
distribution. If a state cannot be written as in Eq. (1) it 
is entangled. 

For a given experimental ^, it is not straightforward to 
prove that it is separable; finding an explicit decomposi- 
tion as in Eq. (1) is a hopeless task. In order to devise 
a simple algorithm for separability testing, we use the 
following two basic facts about separability: 

(i) Consider a separable state ^sep and two other states 
^1 and ^2 such that 



^1 =W2 + (1-P)^sep 



(2) 



for some p G [0; 1]. In this situation, if we can prove that 
g2 is separable, then ^1 must be separable, too. This fol- 
lows directly from the definition in Eq. (1), as this defini- 
tion implies that the convex combination of two separable 
states is separable. Note that if ^2 is entangled this does 
not mean that ^1 is entangled. 

(ii) If a state g is close to the maximally mixed state 
^m = 'Wm^'> ^^^^ i^ is separable. Of course, a pre- 
cise statement of this kind requires a specification of a 
distance. For this, different distances have been inves- 
tigated and bounds on the distance to g^ have been 
obtained (Braunstein et al, 1999; Gurvits and Barnum, 
2002, 2005; Hildebrand, 2006; Kendon et aL, 2002; de Vi- 
cente, 2007; Zyczkowski et al.^ 1998). 

For our purposes, we will use Tr{g^) as a measure of 
the mixedness and use the fact that if 



Trio") < 



1 



NM -I 



(3) 



then g is separable (Gurvits and Barnum, 2002). The 
maximally mixed state has Tr{g^ — 1/NM and states 
with small Tr{g^) are very mixed and close to the maxi- 
mally mixed state. 

Starting from these two facts, the basic idea of our 
algorithm is to find a sequence ^i,^2,^3,--- such that 



Eq. (2) holds for all i,z + 1 and that Tr{g1) > Tr{g1^-^) 
holds. This sequence will lead to a minimization of 
Tr{Ql) and finally, Eq. (3) can be applied to prove sepa- 
rability of Qi and hence of ^i . 

In the simplest case, our algorithm can be sketched as 
follows: 

1. Take the given experimental data ^exp as Qi with 
i = 1. 



2. Consider the optimization problem 



, max 1(01^^10)1 

|0) = |a)|6) 



(4) 



and find a product state \(j)i) = \ai)\hi) which has a 
high overlap with g^. 

Here, the aim is only to find a state with high over- 
lap. In order to make the algorithm work one does 
not need a certified optimal solution of the maxi- 
mization in Eq. (4) (see also below). 



3. Find an Sj with < £,• < Sr, 



such that 



Qi+i '= (l + ^i)^i -^i|0i)(0i| 



(5) 



has no negative eigenvalues and that furthermore 
Tr{gf) > Tr{Qlj^^) holds. Is is natural to choose 
Si such that Tr(^?^]^) is minimal; this is, however, 
not mandatory. 

The point here is that if |0i) has a high overlap with 
Qi^ then it will also have a high overlap with the 
eigenvector corresponding to the maximal eigen- 
value Amax(^i) of Qi. Then, the construction of ^^+1 
is such that usually Amax(^i+i) < Amax(^i) and, 
due to the normalization \min{Qi+i) > Amin(^i) 
holds. Hence, ^i+i will be closer to the maximally 
mixed state than Qi. 

4. Check, whether ^^+i fulfills Eq. (3). If this is the 
case, then ^^+1 is separable and due to Eq. (2) also 
Qi and finally ^exp are separable. Then, the algo- 
rithm can terminate. 

5. If QiJ^i does not fulfil Eq. (3) start again with step 
2 and i ^^ i ^ 1 and iterate further until Eq. (3) 
holds for some i. 

This algorithm deserves some comments: 

(i) First, it is of course not guaranteed that for a sep- 
arable input state the procedure will terminate at some 
point. Consequently, if it does not terminate after many 
steps, one cannot conclude that the state is entangled. It 
is our only claim that the algorithm outlined above is a 
powerful tool in practice. 

In Ref. (Navascues et ai, 2009) an algorithm using 
semidefinite programming has been presented which can 
prove for any non-entangled state the separability after 
a finite number of steps. However, the number of steps 



required is not known in advance, moreover, the applica- 
bility of this technique to larger Hilbert spaces or mul- 
tiparticle problems is not clear. Further algorithms for 
separability testing have been proposed in Refs. (Hulpke 
and BruB, 2005; Spedalieri, 2007); there, however, the 
practical implementation is still missing. 

(ii) Second, note that the algorithm does not require a 
certified solution of any non-trivial optimization problem 
and that it is very robust against imperfections: Even if 
the computation of |0^) or Si is not optimal, this does not 
affect the conclusion that ^exp is separable if Eq. (3) holds 
at some point. Also, the £max is introduced for practical 
purposes, chosing a small Smax makes the convergence of 
the algorithm better in practice. 

(iii) Third, for doing the optimization in Eq. (4) note 
that the optimal |0) = \a)\b) fulfills that \a) is the eigen- 
vector corresponding to the maximal eigenvalue of Xa = 
TrB^Qi^ |^)(^|) and \h) is similarly vector correspond- 
ing to the maximal eigenvalue of X^ = Tr A{Qi\o) {ci\®^) • 
This can be used to tackle the maximization iteratively: 
Starting from a random \a) one computes the optimal \h) 
via Xb^ then with this \h) the optimal \a')^ then again 
the optimal \h') etc. In practice, this converges quickly 
towards the desired solution. 

(iv) Then, one may also check during the iteration 
whether Qi violates some of the usual entanglement cri- 
teria, e.g. the criterion of the positivity of the partial 
transpose. If this is the case, the iteration will never end 
and one can directly stop it. In this case, however, it is 
wrong to conclude that ^exp was entangled. 

(v) Moreover, if one wishes to find an explicit separable 
decomposition of q as in Eq. (1) the algorithm can also 
help: For some of the states close to the maximally states 
explicit decompositions into product vectors are known 
(Braunstein et al.^ 1999). Then, one can write down a 
decomposition consisting of this decomposition and the 
10^) in the iteration. 

The algorithm outlined above can be directly imple- 
mented with few lines of code and performs very well for 
the states generated in the experiment. 



B. Proving separability for the multipartite case 

Let us now demonstrate that the ideas from above can 
also be used to analyze multipartite entanglement. Be- 
fore doing so, we explain some basic notions about the 
entanglement properties of three qubits. 

A pure three-qubit state l?/^) is fully separable, if it is 
of the form l?/^) = |a)|/3)|7) and it is biseparable, if there 
is a grouping of the three parties A, 5, and C in two 
partitions (e.g. AB\C)^ such that it is separable with 
respect to this partition (e.g. I?/^) = \x)ab\v)c)' Other- 
wise, it is genuine tripartite entangled. For the special 
case of three qubits there are further two different classes 
of genuine tripartite entanglement: the GHZ class (repre- 
sented by the GHZ state \GHZ2) = (|000) + |111))/V^) 
and the W class (represented by the W state IVF3) = 



(|001) + |010) + |100))/v^). These states represent dif- 
ferent classes of entanglement in the sense that a single 
copy of a W-type state cannot be converted locally in a 
GHZ-type state (and vice versa), even if the conversion is 
allowed to work only with a small probability (Diir et a/., 
2000). 

As in the bipartite case, one can extend this classifica- 
tion to mixed states via convex combinations: A mixed 
state is fully separable, if it can be written as a mixture of 
fully separable pure states. It is biseparable, if it can be 
written as a mixture of biseparable (and fully separable) 
states. Otherwise, it is genuine multipartite entangled. 
Here, it is important to note that the mixture of bisep- 
arable states may contain pure biseparable states which 
are biseperable with respect to different partitions. Fur- 
thermore, one can define classes of mixed W states and 
GHZ states (Acm et al, 2001). 

For a generalization of the separability algorithm it is 
important that similar results as Eq. (3) exist also for 
multipartite systems (Braunstein et al.^ 1999; Gurvits 
and Barnum, 2005; Hildebrand, 2006; Kendon et al.^ 
2002). For instance, it has been shown in Ref. (Hilde- 
brand, 2006) that if a three-qubit state fulfills 
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Trig") <^^ 0.1397, 



(6) 



then Q is fully separable. More generally, an N-qubit 
state with TV > 3 is fully separable, if 



Tr{Q^) < 



1 



}N 



with a^ 



iN 
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3iv- 



1 



(7) 



hold. This follows also from the results of Ref. (Hilde- 
brand, 2006), where a bound on the radius of the separa- 
ble ball of unnormalized density matrices has been given, 
a rescaling of it delivers Eq. (7). 

Given these facts one can now directly write down 
algorithms to prove full separability or biseparability. 
Finding the fully separable (or biseparable) state with 
the highest overlap (see Eq. (4)) can be done as before, 
by starting from a random fully separable (or bisepera- 
ble) state and then updating it iteratively (see point (iii) 
above). Similarly, one can also write an algorithm which 
can prove that a three-qubit state belongs to the W class, 
since the pure W states can be explicitely parameterized 
(Acm et al, 2000). 

Again, all these algorithms can be easily implemented 
and they have properties similar to the one discussed 
before. 



A. Bipartite distillability 

For the one qubit vs. three qubits partitions (2 x 8), 
it is known that a quantum state in a 2 x A^ system 
is distillable, if it has a negative partial transpose (Diir 
et al.^ 2000). Therefore, distillability in this system can 
directly be checked. 

For the two qubits vs. two qubits partitions (4 x 4) we 
used the mathematical definition of distillability, which 
states that a state g on 1-La ^ ^b is distillable, iff for 
some k we can find four states 

lei) e Ha^-^Ua for zG {1,2} 

k times 

I/,) e Kb^.^^Kb fori e {1,2} (8) 

k times 

such that for \iIj) = a|ei)|/i) + /3| 62)1/2) the estimate 

(^|(^^"f'|^)<0. (9) 

holds (Horodecki et al.^ 1998). For /c = 1 we searched 
numerically for the desired 1?/^), and if this can be found, 
the state must be distillable. 



B. Entanglement superactivation 

Entanglement superactivation is a quantum informa- 
tion processing task, which can work also for the Smolin 
state as a bound entangled state. Consider two copies 
of a four-qubit state g which is distributed among five 
parties as 



^total 



ABCD 



ABCE 



(10) 



that is, the parties A, 5, and C hold two qubits each, 
while the parties D and E have only one. In this situa- 
tion, it was shown that if ^ is a Smolin state, then the five 
parties can create a Bell state between D and E by local 
operations and classical communication only (Shor et a/., 
2003). The protocol uses a sequence of teleportations. 

For the experimental data, we have first applied ap- 
propriate local unitary rotations, then this protocol and 
finally checked with the PPT criterion whether the re- 
sulting state between D and E is entangled. We also 
tested this for arbitrary permutations of the four qubits, 
since the experimental data are, in contrast to the ideal 
Smolin state, not permutationally invariant. 



II. PROVING DISTILLABILITY OF QUANTUM STATES 

Let us now explain the criteria that were used for prov- 
ing distillability of the experimentally generated quan- 
tum states. 



III. CHSH-TYPE BELL INEQUALITY 

We consider the scenario where each of the 4 parties, 
labeled j {j = 1, 2, 3, 4), can choose between two observ- 
ables {of), kj = 1,2. The CHSH-type Beh inequalities 
considered here have the form (Augusiak and Horodecki, 



2006; Werner and Wolf, 2001; Zukowski and Brukner, 
2002): 

\E{1, 1, 1, 1)+E{1, 1, 1, 2)+E{2, 2, 2, 1)-E{2, 2, 2, 2)| < 2 , 

(11) 
where the correlation function E for the measured quan- 
tum state p is an average calculated as follows: 

E{h, fc2, ks, k^Kg) = Tr [pOi\^ Olf ® Olf ® 0["/ . 

(12) 
For the measured states we maximize the left-hand side 
of the inequality in Eq. 11, also known as Bell parameter, 
by varying the observable directions. For the state with 
7 = 0, this Bell parameter is 2.21(2), while for the state 
with 7 = 0.06, the parameter is already 1.47(3). 



IV. BISEPARABILITY OF THE STATES WITH 

7 = 0,0.06,0.12, OR 0.18. 

Applying the separability algorithm to the measured 
states for these values of 7 does not lead to a proof that 
these states are biseparable. This is due to the fact that 
the experimental states have some eigenvalues which are 
practically zero, making the third step of the separability 
algorithm (Eq. 5) difficult. However, if one takes the state 
with 7 = and adds 2% of white noise 



0.98^(7 = 0) +0.02 



1 
16 



(13) 



the algorithm proves that this state ^noise is bisepara- 
ble. It should be noted that the statistical fluctuations 
of the diagonal elements of ^(7 = 0) are larger than the 
amount of noise added. Moreover, if one considers the 
state consisting of the mixture of all Monte Carlo sample 
states, the algorithm proves that this state is bisepara- 
ble. This gives strong evidence for the biseparability of 
the experimentally generated states. 



V. MODEL OF IMPERFECTIONS REVEALS THE 
OBSERVED DYNAMICS 



State preparation imperfections: 
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FIG. SI Negativity for the 2:2 bipartition (left axis) and 1:3 
bipartition (right axis) as a function of partial decoherence for 
the shown imperfections (see text). The imperfections have 
no effect on the dynamics of Ni:3. The negativity is the same 
for all permutations of each bipartition due to the symmetry 
of the initial state. 



VI. QUANTUM STATE RECONSTRUCTION AND 
STATISTICS 

Physical states were tomographically reconstructed via 
maximum likelihood estimation (MLE) (Jezek et a/., 
2003) . On the four qubits, we performed 3x3x3x3 = 81 
measurements by iterating the measurement of the ex- 
pectation values of ax^ (Jy^ and (Tz on each qubit. For 
each state, the number of copies used per state are shown 
in Table L Tomographic data were acquired in less than 
5.5 hours for each state. Full tomographic sets of 200-250 
copies per measurement were acquired iteratively to keep 
track of potential drifts in the experimental setup. 

Error analysis was calculated via Monte Carlo simu- 
lations over the multinomially distributed measurement 
outcomes of the state tomography. For each state, 200 
Monte Carlo samples were generated and reconstructed 
via MLE. 



Several imperfections affect the experimental setup, 
but to briefly demonstrate that the experimentally ob- 
served entanglement dynamics agrees with straightfor- 
ward calculations, here we choose to only assume im- 
perfect M0lmer-S0rensen operation times (angle in 
MS{0^(j))^ see report). We consider longer operation 
times by a small fraction, si for the 3-out-of-4-entangling 
operation, MS{{1 + £:i)7r/4,0), and £2 for the GHZ- 
entangling operation, MS{{1 + £:2)7r/2,0). Using the 
preparation sequence shown in Figure 1 of our report 
and assuming other operations are perfectly realized, we 
calculate an expected initial state. Upon decohering this 
imperfect initial state in two steps, as described in the 
report, we observe the dynamics of entanglement resem- 
bling our measurements, as illustrated here in Figure 1. 
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TABLE I Negativity and 


smallest eigenvalue of the partial transpose of measured states for every bipartition. 








2:2 & 1:3 entangled 


Bound Entangled 


Fully Separable 


single"^ 


7 = 


0.00 


0.06 0.12 0.18 


0.24 


0.32 


0.35 


0.47 


0.60 


0.43 




ABiCD 


0.033(5) 


0.046(6) 0.017(5) 0.009(5) 


0.000(3) 

















N2:2 


ACiBD 


0.041(5) 


0.043(6) 0.012(5) 0.007(5) 


0.000(3) 



















AD:BC 


0.044(5) 


0.047(7) 0.013(5) 0.006(5) 


0.000(3) 



















AB:CD 


-0.011(2) 


-0.016(2) -0.008(2) -0.005(2) 


0.005(2) 


0.020(2) 


0.018(3) 


0.028(3) 


0.043(2) 


0.013(2) 


min(eig(p^2:2)) 


AC:BD 


-0.012(2) 


-0.016(2) -0.006(2) -0.004(2) 


0.004(3) 


0.019(2) 


0.021(2) 


0.029(3) 


0.045(2) 


0.011(2) 




ADiBC 


-0.013(1) 


-0.015(2) -0.007(2) -0.003(2) 


0.005(3) 


0.022(2) 


0.016(3) 


0.031(2) 


0.043(2) 


0.013(2) 




A:BCD 


0.715(8) 


0.507(14) 0.438(14) 0.335(13) 


0.199(11) 


0.035(7) 


0.021(10) 








0.038(7) 


iVl:3 


BiACD 


0.715(8) 


0.509(14) 0.440(14) 0.337(12) 


0.195(11) 


0.032(8) 


0.022(10) 








0.039(7) 




CiABD 


0.715(8) 


0.510(14) 0.439(14) 0.337(12) 


0.197(11) 


0.038(8) 


0.015(9) 








0.042(7) 




D:ABC 


0.716(8) 


0.509(14) 0.438(14) 0.337(12) 


0.208(11) 


0.045(7) 


0.028(10) 








0.045(7) 




A:BCD 


-0.102(2) 


-0.088(2) -0.077(3) -0.064(3) 


-0.045(3) 


-0.015(2) 


-0.008(3) 


0.019(3) 


0.042(2) 


-0.011(2) 


min(eig(p^2:2)) 


BiACD 


-0.104(2) 


-0.088(3) -0.079(3) -0.063(3) 


-0.040(3) 


-0.014(2) 


-0.008(3) 


0.020(3) 


0.042(2) 


-0.015(2) 




CiABD 


-0.104(1) 


-0.087(3) -0.075(3) -0.062(3) 


-0.042(3) 


-0.012(2) 


-0.007(3) 


0.019(3) 


0.042(2) 


-0.014(2) 




D:ABC 


-0.105(2) 


-0.089(2) -0.079(3) -0.064(3) 


-0.045(3) 


-0.016(2) 


-0.009(3) 


0.018(3) 


0.043(2) 


-0.014(2) 


(state copies)/ (81 meas.) 


5000 


3000 3000 3250 


3750 


6000 


3100 


3000 


4000 


6000 



^Bound-entangled state decohered in a single step, see report. 



TABLE II Summary of Monte Carlo samples satisfying the tested properties. 



Decoherence 7 = 





0.06 


0.12 


0.18 


0.24 


0.32 


0.35 


0.47 


0.60 


single 


2:2 entangled 


200 


200 


200 


200 


41- 














1 


2:2 distillable 


200 


200 


194 
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2^ 

















1:3 entangled 


200 


200 


200 


200 


200 


200 


200 








200 


1:3 distillable 


200 


200 


200 


200 


200 


200 


200 








200 


superactivatable 


200 


200 


200 























biseparable 


200" 


200" 


200" 


200" 


200 


200 


200 


200 


200 


200 


fully separable 


























200 






^These 41 states had a non-positive partial transpose for all per- 
mutations of the 2:2 bipartition. 

^44 states had a positive partial transpose on every permutation 
of the 2:2 bipartition (undistillable) and were also 2:2 separable. 

^Strong evidence, discussed in text. 



TABLE III Absolute value of measured density matrices. 

|(p(7)|for7 = |(p(7)|for7 = 0.06 
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TABLE IV Absolute value of density matrices ^(7) calculated by applying local unitaries to the measured states ^(7) which 
maximize the amplitudes of their real part. This transformation facilitates visualizing the Smolin components of the states. 

I (£.(7) I for 7 = 1(^(7)1 for 7 = 0.06 | (£.(7)! for 7 = 0.12 




1(^(7)1 for 7 = 0.18 




|(e(7)|for7 = 0.24 



|(e(7)|for7 = 0.32 




1(^(7)1 for 7 = 0.35 



|(e(7)|for7 = 0.47 



|(e(7)| for 7 = 0.60 




s'^^ #.^ 




0000 



I(f?(7)|for7' = 0.43 




